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Abstract: Perpendicular electron dynamics and the associated collisions are discussed 
in relation to the collisional drift wave instability. In addition, the limit of small parallel 
wave numbers of this instability is studied and it is shown to yield a reduced wave fre- 
quency. It is also shown that in this case the growth rate in fact decreases for smaller 
parallel wave numbers, instead of growing proportional to l/A;^. As a result, the growth 
rate appears to be angle dependent and to reach a maximum for some specific direction 
of propagation. The explanation for this strange behavior is given. A similar analysis is 
performed for acoustic perturbations in plasmas with unmagnetized ions and magnetized 
electrons, in the presence of a density gradient. 



PACS Numbers: 52.35.Kt, 52.30.Ex, 52.20.Fs, 52.35.Fp 



I. INTRODUCTION AND MODEL 



Typically, collisions in plasmas are responsible for the damping of plasma modes. 
However, this does not hold for the drift wave which grows in the presence of electron 
collisions. The growth appears because, in the presence of collisions, the perturbed wave 
potential lags behind the perturbed density. This, together with the free energy stored 
in the density gradient, makes the mode growing. In the case when the ions are not 
magnetized (e.g., perturbations with frequencies above the ion gyrofrequency Qi, or /and 
when the ion collision frequency is above fli), the ion motion will be of the sound-type, 
while the motion of electrons can still be described within the drift wave limits. Hence, 
acoustic-type perturbations in such a case may also become unstable due to the same 
reasons, ^'^ with the instability taking place provided that the acoustic wave frequency is 
&e/ou;the electron diamagnetic frequency cu^e — v^eky, v*e = ~{''^Te/^e)^z x Vno/rio, v^^ — 
KTg/rrie. Physically, the instability is due to the fact that collisions prevent the electrons 
from moving freely and from shielding the electrostatic perturbations that involve the ion 
component. 

In most cases, these electron collision effects are included from the electron parallel 
momentum. This is because of the small inertia, which results in a dominant electron 
motion in the parallel direction, as compared to the perpendicular gyromotion-type dy- 
namics. Yet, a clear distinction about when such a model is justified is rarely seen in the 
literature. In our recent work^ it is shown that the model can be used provided that 

k'M/ikye) » 1- (1) 

Here, kz and ky are the wave number components in the direction parallel and perpendic- 
ular to the direction of the magnetic field vector Bq = BqCz, respectively, and He denotes 
the electron gyrofrequency. The electron collision frequency Ve may include collisions with 
both ions and neutrals. Numerous parameters contribute to the condition (1), and conse- 
quently it may not always be satisfied, so that in some cases it may become necessary to 
include collisions in the perpendicular electron dynamics as well. An obvious consequence 
of these coUisions is that the usual drift motion of electrons (perpendicular to the force 
vector) will be modified, and an additional component of motion in the direction of force 
vector will appear as well. 

In a weakly ionized plasma, the electrons may predominantly collide with neutrals. 
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Usually, these collisions are described assuming a static and heavy neutral background so 
that the dynamics of the neutrals is ignored. However, the interaction between the two 
fluids happens due to friction, and the momentum conservation imposes that the neutral 
dynamics must be taken into account in spite of a big difference of mass (per unit volume) 
of the two fluids. In a recent study^ the differences between the two models are discussed 
in detail. 

In the present work, we consistently take the dynamics of the target particles (neutrals 
in the present case) into account, and we discuss the domain in which the condition (1) 
is not satisfied. 

In the case when the temperature of the heavy plasma constituents (ions and neutrals) 
is much below the electron temperature, we omit their thermal effects. The electron 
momentum equation is 

= eneV0 - eUeVe X B - KTeVUe - nieneiyeniVe " Vn), (2) 

and for the neutrals we use 

Vn = -l^ne{Vn- Ve). (3) 

The momentum conservation implies that = rneTiei'en/ {'iT^rJ^n) ■ 

II. ON THE EQUILIBRIUM 

It is seen that, assuming a stationary equilibrium and also {vnQ-V)vnQ — 0, {veQ-V)veQ — 
0, from (3) we have v^o — '^en- At the same time, Eq. (2) yields Veo — v*e- Hence, due to 
friction both fluids move together. In most cases under laboratory conditions such a state 
is not expected to take place because of the long time that is needed to reach it. In order 
to get a better feeling about the necessary time scales, we may set a constant velocity for 
the electrons Veo (assuming that the parameters determining the velocity are controlled 
externally). Using Eq. (3) we then obtain the neutral velocity: fno/^eo = 1 ~ cxp(— z/„et). 
Now, setting for example Ueo = 10^^ m~^, and assuming a hydrogen plasma in hydrogen 
gas with Uno = 10^^ m~^, Te = 1 eV, we have — 1-09 • 10^ Hz, z/„e = 0.6 Hz, and we 
find that Vno/veo — 1 within about 8 s. Here i^en = CTennnoVre, and^ aen{Te) — 2.6 • 10~^^ 
m^. 



meUe 



dVe 

dt 
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Similarly, setting rieo = 10^''' m"^, n„o = 10^^ m"^, Tg = 5 eV, we have i^en = 1-07 • 10^ 
Hz, — 5.8 Hz, and Vno/veo — 1 within about 0.8 s. Here^ (7e„(Te) ~ 1.14 • 10~^^ m^. 

Numbers used here are for demonstration only. Yet, similar parameters may be seen 
in various studies dealing with coUisional drift-waves, like in the experimental work from 
Ref. 9 (with = 2 eV, < 0.1 eV, Ueo = 10^^ m-^ and i^en = 2 • 10^ Hz). A similar 
density (Ueo = 7 ■ 10^^ m^^) but a 5 times lower electron temperature is used in the 
experiment in Ref. 10, while in the experiment in Ref. 11 the density and the temperature 
arc also of the same order as in our examples, i.e., rieo = 5 ■ 10^"'^ m^^, Tg ~ 1 eV. In 
another experiments^ the parameters are Ueo = 2 • 10^^ m~^, Te — 2 eV. 

Although parameters may have very different values in various experiments, the num- 
ber used in the examples above are frequently seen in the laboratory plasmas. Although 
the ionization degree in the two examples given above is not so small, the electron-neutral 
collisions are larger than electron-ion collisions. For the given parameters, the two cases 
yield (in kHz) i/^i ~ 304, 305 < i/gn, respectively. Here ly^i remains almost the same be- 
cause both the density and the temperature are increased. We used^^'^^ the expression 
for the collision frequency between any two charged species b and a 



Here, L^a — log{Xd/bo) is the Coulomb logarithm, = Xdb^da/ {Xdb + ^daY^"^ is the plasma 
Debye radius, and bo — [\qbqa\ /[12TreoK{Th + Ta)] is the impact parameter for Coulomb 
collisions. 

Hence, it is correct to say that in most cases we have a 'time evolving background 
plasma'. Yet, this time dependence is negligible as long as the period of the perturbations 
is much shorter than the equilibrium evolution time (that is of the order of l/z/„e). In 
fact, this is not a rough approximation bearing in mind the other plasma parameters 
and the corresponding evolution. For example, the assumed equilibrium density gradient 
is also time dependent because of the electron collisions, and those take place on much 
shorter time scales. In other words, fast electron collisions will tend to destroy the density 
gradient, and this may happen on short time scales, unless the density gradient is somehow 
kept externally (examples of that kind can be found in Refs. 15, 16). This can be seen by 
taking the second set of the above given parameters and assuming Bq — 0.01 T, and the 
plasma density gradient length as L„ = [{dueo/ dx) /neo]~^ = 0.1 m. The perpendicular 
electron diffusion velocity is then given by Dg/Ln = 3 m/s, where Dg = z/en^re/^e- Hence, 
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electron diffusion perpendicular to the magnetic field vector will tend to destroy the given 
density gradient in about 0.03 s and all this is entirely due to collisions. Regarding the 
neutrals, we stress that there may also be a continuous diffusion of neutrals from nearby 
regions, and those new particles have no time to be accelerated by friction within the 
assumed short wave period. As a result, in most cases one may perform a wave analysis 
neglecting the equilibrium movement of the neutrals. 



III. DERIVATIONS AND RESULTS 
A. Drift wave 

From Eq. (2), in the limit when the wave phase velocity and the perturbed electron 
velocity are much below the electron thermal velocity (equivalent to the usual massless 
electron approximation in the literature), we obtain the electron perpendicular velocity 



1 

Ve± 



1 + u^^ayni 



— X V_L0 + — r ^-7^ TT^e, X 



• (5) 



Bo He Bo Ql He rig 

Here, a = uj/{uj + iVne)- In view of the previous discussion on equilibrium, the perturbed 
neutral velocity used here is Vn = WneVe/i^ + ii^ne)- The electron parallel momentum 
yields the parallel velocity 



Using Eqs. (5, 6), and neglecting |z/g„a^/f2g| in comparison to 1 in the denominator in 
Eq. (5) (the limit of magnetized electrons), from the electron continuity equation 

duei 



we obtain 



Here 



+ V±(neiT±e) + ^zirieOVezl) = 0, 

Ueo uj + iDp + iD^{uj'^ + ^'^J/(a;2 - ivne^) i^T^ ' 



(7) 



Dp = VenOlklpl, = klvljPen: Pe = Vre/^e- 

The term Dp describes the effects of electron coUisions in the perpendicular direction. 
Note that neglecting the neutral dynamics here is equivalent to setting — 0, yielding 
a = 1, and this corresponds to the model in Refs. 1, 2. In the case Dp = 0, Eq. (7) 
becomes the same as the corresponding equation from Ref. 6. It is seen also that the 
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assumption \Dz/Dp\ ^ 1 yields the condition (1). In the present work we are interested 
in the domain when this hmit is either reversed or the ratio is of order of 1. 

For cold coUisionless ions we use the same model as in Ref. 6, i.e., the ions are only 
subject to the electromagnetic force, so that the momentum equation for magnetized ions 
is 

miUi 



—eniV(p + eniVi x B. (8) 
Using this, from the ion continuity we obtain 



+ P. = c.M, cl^KTjm,. (9) 

The dispersion equation within the quasi-neutrality limit reads 

The term with the sound speed is due to the ion parallel response and it is known to 
introduce the threshold for the instabihty.^ 

For negligible collisious, the right-hand side in Eq. (10) reduces to 1, and we have the 
simple expression for the drift wave in an ideal plasma 



[l + klpl)-uj.eOJ-klcl = 0. (11) 



It has two solutions, one positive coupled drift-acoustic mode, and one negative acoustic 
mode. 

The full Eq. (10) is solved numerically for several values of the neutral density. In 
Fig. 1 the frequency (normalized to cu^e) of the positive drift-acoustic solution is given 
in terms of the parallel wave-number kz, for a fixed ky = 10'^ l/m, and for Uno = 10^^ 
neo = 5 ■ 10^6 m-^ Te = 11600 K, Bq = 0.1 T, and L„ = 0.1 m. This particular 
set of parameters can easily be achieved in laboratory conditions, and it is taken in 
order to have perfectly satisfied all conditions used in the model. This is seen from the 
following, cu^e — 10^ Hz, and the mode frequency is below cu^e or of the same order, 
so that for the given parameters both electrons and ions are magnetized, Jl^ ~ 10^ Hz, 
Jle = 2 • lO^'^ Hz. The electron coUision frequency with ions is about 3 orders of 
magnitude below — 10^ Hz, and it is therefore negligible. The plasma P is about 10"^ 
and the electrostatic limit is justified. Large ky impfies a reasonably well satisfied local 
analysis condition because Ln/Xy ~ 16. In the range of frequencies and kz given in Fig. 1, 



for D=0 
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p 
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Figure 1: The positive solutions of Eq. (10) without (dashed hne) and with (fuU hnc) the 
perpendicular electron collisions included. Dotted line is the positive solution of Eq. (11). 
Frequency is normalized to a;*e- 

the ratio {uj/kz)/vTe changes between 0.12 (for kz = 1) and 0.02 (for kz = 14), so that 
neglecting the left-hand side of the electron momentum equation is justified. Observe 
also that the electron mean free path is Vre/j^en — 0.004 m, while the shortest parallel 
wavelength is about 0.45 m, so the condition for the fluid model for electrons is well 
satisfied. 

The dotted line in Fig. 1 is the solution which follows from the ideal case (11) (i.e., for 
the right-hand side in Eq. (10) equal to unity). In the hmit of small kz this solution goes 
to 0.49 a;*e because here kyPs = 1.02. The dashed line is the numerical solution of Eq. (10) 
in the case when the electron collisions in the perpendicular direction are simply ignored 
{Dp = 0). The full line is the solution including Dp. Observe the significant difference in 
the limit of small kz between the ideal mode frequency and the collisional modes. The 
effect of keeping Dp here is practically to restore (at least quantitatively) the drift-like 
behavior of the mode. 

The growth-rates associated with the frequency from Fig. 1 are presented in Fig. 2, 
showing an angle dependence and having a maximum at kz around 3 (for the case Dp — 0). 
We note that a similar angle dependence (though due to ion coUisions) was obtained 
recently in the problem of current driven acoustic modes. ^'^^ In the present case, the 
origin of the maximum, and of the reduced frequency for small kz is different and can be 
understood from the following. If Dp is completely omitted, the approximative growth 
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Figure 2: Growth rate (normalized to ou^g) associated with the mode from Fig. 1. The 
dash-dot hne refers to the growth rate which follows from Eq. (12), i.e., in the limits 
D„ = Q and bl < D,). 



rate is given by^ 



CO, 



2 i^fky, - kici 



(12) 



Here uir,i denote the real and imaginary parts of the frequency, respectively. This analytic 
expression is obtained and valid within the approximation 



\(jJr\ < D, 



(13) 



Hence, in this limit, the growth rate appears to be proportional to l/D^ ~ l/Zc^ and it 
strongly grows for small k^. This is presented by the dash-dotted line in Fig. 2. 
The condition (13) can be written as 



\uJr\/kz 



kyVrj 



(14) 



The right-hand side is the ratio of the mean free path and the parallel wavelength, and 
it must be below 1 to have a proper fluid theory. The left-hand side (the ratio of the 
parallel wave phase speed and the electron thermal speed) must also be below 1 within the 
model used here (neglected left-hand side of the electron momentum equation). However, 
although both sides separately have some physical meaning and both are below 1, their 
mutual ratio is arbitrary and, therefore, the conditions (13,14) do not have to be satisfied. 
In fact, in our present case, for the parameters used above and in the limit of small k^, 
the frequency \ujr\ becomes of the order of or even larger. This reverses the direction 
of the dash-dotted line in Fig. 2 for small enough fc^, and we obtain the dashed line (the 
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Figure 3: Dp — Uenkypl normalized to in terms of kz corresponding to Figs. 1, 2. 

case Dp — 0). Hence, for small kz the growth rate takes a completely different shape if 
the condition (1) is satisfied, and if the dispersion equation is solved for parameters when 
the condition (13) does not hold. 

Remark that the condition (13) is usually used for analytical convenience. We have 
shown here that this excludes an important limit in which the growth rate may be dras- 
tically reduced. 

In addition to this, in the limit opposite to the condition (1), the perpendicular electron 
dynamics and perpendicular collisions must be taken into account. As a result the growth 
rate is additionally modified and it is given by the full line in Fig. 2. This effect appears 
because for /c^ < 6 1/m, we have Dp > D^ as seen from Fig. 3, where Dp is presented 
normalized to local values of Dz{kz). Basically the same reasons are behind the reduced 
frequency in Fig. 1. 

Note also that the growth rate in Fig. 2 becomes negative around kz ^ 10.45 1/m. 
The reason for this may be seen from Eq. (12), where the growth rate changes the sign 
if fc^ > u'^kyp'^/cj.. This is the instability threshold mentioned earlier, due to the parallel 
ion response. The point of the sign change is almost the same for all three cases from 
Fig. 2 because all parameters, except the frequency, are kept constant, and the frequency 
itself for the three cases is almost the same (see Fig. 1). 

Equation (10) is solved also in terms of ky, for a fixed kz — 4: and for the same 
values for the other parameters as above. This yields Dz — 25810 Hz and ps — 10~^ m. 
The result for the frequency and the growth rate (both in units of cu^e) is presented in 

— * 

Figs. 4, 5. In the given range of k the ratio Dp/Dz has values in the range from 0.2 (for 
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Figure 4: The normalized frequency as positive solution of Eq. (10) in terms of ky. 
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ky — 300 1/m) to 15 (for ky — 2500). The growth rate changes the sign at kyPs ^ 0.65. 
Setting Dp — provides the dashed hnes in Figs. 4 and 5. The growth rate is lower in 
the presence of Dp. Physically, in the situation with Dp ^ 0, hght electrons have more 
possibility to shield electrostatic perturbations by moving now also in the perpendicular 
direction (due to the collisions). As a result, the growth rate is reduced. 

B. Acoustic wave 

In the case of high wave frequencies 

> 1^1 > ^i, (15) 

the ions are unmagnetized and their response is the same as in an ion acoustic mode, 
while the electron dynamics is the same as in the previous text. The ion continuity and 
momentum equation (8) without the Lorentz force instead of Eq. (9) yield nn/riiQ — 
(/c^c^/cc;^)e0i/(fi;Te), k'^ = ky + k^. The dispersion equation now reads 

cu^ u + iDp + iDz{u!'^ + 1/2^) /(cj2 _ iu^^ui) 

Equation (16) is solved for modified parameters in order to have the conditions (15) 
satisfied. Hence, we take heavier (argon) ions rrii — AOrUp, rUn — rrii, and Tg — 5 eV, 
Uno — 10^^ m~^, Ueo — — 5 ■ 10^^ m"^, Bq — 0.01 T, L„ = 0.1 m. For these parameters 
we have^^ = 8.7 • 10~^° m^. We set ky ^ 7 ■ 10^ 1/m and solve Eq. (16) in terms 
of kz- The results are given in Figs. 6 and 7. It is seen that the perpendicular electron 
collisions drastically destabilize the mode. In the limit of small kz, the growth rate in 
Fig. 7 is about 70 times larger. Note that for k^ = 0.3 1/m we have Dp/D^ = 141, while 
for A;^ = 14 1/m this ratio is only 0.06. 

For large k^ values, the growth rate in Fig. 7 decreases with k^ similar to the solutions 
from Rcf. 2. It can also be shown that the growth rate becomes larger when n„o is increased 
(i.e., for larger collision frequency). These features are in agreement with an approximative 
analytical solution that may be derived from Eq. (16) assuming a dominance of the D^ 
term. 

IV. CONCLUSIONS 

To summarize, we have presented some details of the coUisional drift wave instability 
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Figure 6: Ion acoustic frequency (normalized to a;*e) as the positive solution of Eq. (16). 
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Figure 7: The normalized growth rate associated with the acoustic mode from Fig. 6. 
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that are usually overlooked in the literature. In particular, the mode behavior for small 
parallel wave numbers is usually described within the limit given by the condition (13), 
resulting in a strongly increasing growth rate proportional to We have shown that in 
the opposite limit the mode is considerably modified, both the frequency and the growth 
rate are reduced. As a result, the growth rate appears to be angle dependent, reaching a 
maximum value for a certain angle of propagation. Furthermore, the effect of collisions 
on the perpendicular electron dynamics and on the collisional drift wave instability is 
discussed. It is shown to have effect on both the wave frequency and the growth rate, as 
presented in Figs. 1, 2, 4, and 5. 

Similar effects of the perpendicular electron collisions are shown for the acoustic-type 
mode in plasmas with magnetized electrons and unmagnetized ions. In such plasmas, 
the ion electrostatic perturbations result in acoustic modes regardless of the angle of 
propagation. The physics of the mode is described in the classic Ref. 1, and was also 
discussed more recently in Ref. 2. 

The collision frequency has been taken constant in the present work, following standard 
semi-empirical models used in the literature, yet a recent study^^ shows that it can be 
significantly modified depending on ion temperature. This is an issue that requires further 
investigations. In any case we believe that the results presented here should be taken into 
account in studies dealing with dissipative drift-type instabilities, like those in Refs. 20-23. 
A most obvious example where the instability conditions will clearly be modified is the 
recent Ref. 24 dealing with drift dissipative instability in a toroidal device. 

The model and results obtained here may be used also in application to space plasmas. 
A most obvious example where the theory works is the lower solar atmosphere with the 
ionization ratio which may be as low as 10~^ (in photosphere, with typical temperatures 
of around 5 ■ 10^ K), and it changes with the altitude to become of the order of 1 in 
the chromosphere (at the altitude of about 2200 km), and grows towards corona where 
neutrals are absent. Some aspects of the drift waves in solar atmosphere may be found 
in Refs. 25, 26. A similar situation with altitude-dependent parameters suitable for the 
drift wave may be found in the terrestrial ionosphere.^^ Another interesting aspect of the 
application to space problems may be seen also in Ref. 28 dealing with Hall thrusters used 
in spacecrafts propulsion, with dominant collisions with neutrals and a geometry suitable 
for the drift wave analysis. 

We have briefiy discussed the issue of the equilibrium in such a three-component 
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coUisional plasma. One important point that follows from it, and that has been discussed 
briefly, is the possibility of equilibrium macroscopic flows of neutrals in a magnetized 
plasma. Such flows may be generated by diamagnetic drifts of plasma species in an 
inhomogeneous plasma. The diamagnetic drift itself is not a macroscopic plasma flow 
but the effect of gyromotion of plasma species in the presence of a density gradient. 
Yet, as such it may still generate real macroscopic flows of the neutral component of the 
plasma. This is due to friction and may have important implications, in particular for 
space plasmas. Similar effects should appear in the presence of a magnetic field gradient. 

Acknowledgements: The results presented here are obtained in the framework of the 
projects G.0304.07 (FWO-Vlaanderen), C 90205 (Prodex), and GOA/2009-009 (K.U. 
Leuven) . 
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